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SUMMARY 


The maximum likelihood parameter estimation technique was used to extract 
values of stability and control derivatives from flight test data obtained from 
a light, single-engine, low-wing airplane. The flight tests consisted of 9 runs 
in which the stabilator was used to excite the longitudinal motions and 28 runs 
in which the rudder and ailerons were used to excite lateral motions. The vari- 
ous control inputs were initiated from trimmed level flight with a trimmed air- 
speed of about 46 m/sec and an initial altitude of about 600 m. 

The consistency of the derivative estimates as it relates to various inputs 
was investigated to determine the inputs which provide the most information for 
identification. Three criteria were used in this investigation: the ensemble 

variance, the estimated Cramer-Rao lower bound, and parameter correlations. 

There were no significant differences in the consistencies obtained by using the 
various stabilator inputs. On the other hand, for the lateral case, the sequen- 
tial inputs (rudder followed by ailerons or ailerons followed by rudder) gave 
noticeable improvement in parameter value consistency over the rudder or aileron 
inputs individually. Also, some improvement in consistency was noted for both 
aileron and rudder square waves as opposed to sine waves. 

The derivative values obtained by using the maximum likelihood technique 
were compared with values computed by using an empirical approach. They com- 
pared favorably in direct comparison and in comparison of time histories of com- 
puted airplane motions, although the extracted parameters provided the better 
match. 


INTRODUCTION 

The maximum likelihood technique (ref. 1) has been used to estimate values 
of the aerodynamic derivatives for several airplanes ranging from experimental 
high-performance airplanes (refs. 2 and 3) to single-engine light airplanes 
(ref. 4). The importance of proper inputs to the success of the parameter iden- 
tification process has long been recognized (ref. 5), and several input design 
studies have been conducted using modern control theory (refs. 5 to 8). The 
inputs developed in these studies can be difficult to implement in practice, and 
approximations of these optimal inputs are often used. The approach taken in 
this paper to determine the effectiveness of inputs in parameter identification 
was to perform repeatedly several different inputs which were roughly based on 
some of these designed inputs (refs. 6 to 8). The resulting parameter estimates 
for each type of input were compared to indicate how well the inputs provided 
the information required to define each parameter precisely. 

The purpose of this paper is to present the evaluation of the effectiveness 
of different inputs in providing consistent parameter estimates. This was done 
by conducting flight tests using several types of inputs, extracting parameters 


from the flight test records, and using three criteria to determine the relative 
effectiveness of the various inputs relating to parameter estimation. The three 
criteria used for input evaluation in the study were to compare ( 1 ) the vari- 
ances of the estimated parameter values, (2) the estimates of the Cramer-Rao 
lower bound (ref. 9), and (3) the parameter correlations. 

This report first describes the airplane, instrumentation, and flight 
tests. Then the control inputs used in the tests are described and the three 
criteria used for comparing these inputs are discussed briefly. The results of 
these comparisons are presented next. Finally, a comparison of the parameter 
estimates with those computed by using an empirical-theoretical approach is 
presented. 


SYMBOLS 

The aerodynamic parameters are referenced to a system of body axes with 
the origin at the airplane center of gravity, which is located at 21.25 percent 
c, and with orientation of body axes as shown in figure 1. 

ax>ay,az acceleration measured along X, Y, and Z body axis, respectively, 
g units 

b wing span, m 

c wing mean geometric chord, m 

Fx,Fy,Fz force along X, Y, and Z body axis, respectively, N 
g acceleration due to gravity, m/sec^ 

Ix,Iy>Iz moment of inertia about X, Y, and Z body axis, respectively, 
kg-m^ 

Ixz product of inertia, kg-m^ 

i index 

J cost function 

L likelihood function 

It distance from airplane center of gravity to center of pressure of 

horizontal tail, m 

Mx,My,Mz rolling, pitching, and yawing moments, respectively, N-m 
m mass, kg 

N number of data points 

n dimension of system 
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roll rate, rad/sec 

number of parameters to be identified 

pitch rate, rad/sec 

dynamic pressure , N/m^ 

estimate of error covariance matrix 

yaw rate, rad/sec 

wing area , m^ 

standard deviation 

standard deviation of estimated Cramer-Rao lower bound 
thrust , N 

velocity along X, Y, and Z body axis, respectively, m/sec 

velocity component along X, Y, and Z body axis, respectively, at 
angle-of-attack sensor on wing-tip boom, m/sec 

airplane total velocity, m/sec 

body coordinate axes through airplane center of gravity 
body coordinates, m 

X-, y-, and z-coordinates , respectively, of the sensors on wing-tip 
boom relative to airplane center of gravity, m 

vector describing state of airplane 

angle of attack, rad 

parameter vector 

parameter change vector 

angle of sideslip, rad 

left aileron deflection minus right aileron deflection, rad 
stabilator deflection, rad 
rudder deflection, rad 
aileron sine-wave input 
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(5a)sq aileron square-wave input 

(6e)pul stabilator rapid pulse followed by slower decay 

(^e^si stabilator sine-wave input 

(6e)sq stabilator square-wave input 

(6r^si rudder sine-wave input 

(6r)sq rudder square-wave input 

(6a-6j.)si aileron sine-wave input followed by rudder sine-wave input 

(6a“^r^sq aileron square-wave input followed by rudder square-wave input 

(6p-6a)si rudder sine-wave input followed by aileron sine-wave input 

(6r«-5a^sq rudder square-wave input followed by aileron square-wave input 

e angle between thrust axis and airplane X body axis, positive for 

thrust up , rad 

ri measurement noise vector (n x 1 ) 

0 pitch angle, rad 

p air density, kg/m3 

an m.R ensemble mean estimated Cramer-Rao lower bound 

(j) roll angle, rad 

rolling-moment coefficient, My/gSb 
Cjjj pitching-moment coefficient, My/dSc 

Cn yawing-moment coefficient, M^/qSb 

Cx thrust coefficient, T/qS (used in some publications as Tq’) 

Cx axial-force coefficient, Fx/qS 

Cy side- force coefficient, Fy/qS 



normal-force coefficient, Fx/qS 
9Cx ^ 9Cx 

3 Pb a £b 

2V 2V 
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Subscripts : 
c computed 

k index 

m measured 

o coefficient at trimmed conditions 

t trimmed conditions 

Superscripts: 

T transpose matrix 

M measured quantity 

o nominal evaluation 

A dot over a symbol signifies a derivative with respect to time. 
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DESCRIPTION OF AIRPLANE, INSTRUMENTATION, FLIGHT TESTS, 

AND DATA REDUCTION 


The subject airplane for this study was a four-place, low-wing, single- 
engine airplane shown in figure 2. Its pertinent geometric details are given 
in table I. The movable control surfaces included the stabilator, rudder, aile- 
rons, and flaps; and the airplane was instrumented to record control-surface 
movements and airplane responses to these movements. 

The variables recorded for this study were obtained from an onboard instru- 
mentation package. The range of each instrument used to record the pertinent 
variables is given in table II. The accelerations, angular rates, and angular 
attitudes were recorded continuously, whereas the airspeed, angle of attack, 
angle of sideslip, and control-surface positions were recorded sequentially by 
use of a commutator which sampled each variable 20 times per second. Airspeed, 
angle of attack, and angle of sideslip were measured on a boom that was located 
near the wing tip and extended 3/4 c ahead of the wing leading edge. (See 
fig. 2.) 

The flight tests were divided into two groups: those in which the longitu- 

dinal modes were excited and those in which the lateral modes were excited. No 
flaps were used during the flight tests and all tests were initiated from trimmed 
level flight. The test procedure was as follows: The pilot turned on the data 

recording equipment, made specific control inputs, allowed airplane responses 
to settle out, and then turned off the recording equipment. This sequence con- 
stituted a data run. The tests were conducted in smooth air to minimize process 
noise from gusts. The throttle was held fixed during the runs to minimize any 
changes in thrust so that the only disturbing force on the airplane was due to 
the control input. Since angle-of-attack oscillations during the tests were 
generally less than 7° peak to peak, these variations were considered minimal, 
although the reader should note possible thrust change effects on and 

as stated in appendix A. 

A total of 37 flight test runs were made, 9 runs concentrating on longitu- 
dinal dynamics and the remainder concentrating on lateral dynamics. The length 
of the runs ranged from about 20 to 50 sec. Since the- object of the flight test 
was to evaluate the relative effectiveness of various inputs for parameter esti- 
mation, all test runs were initiated from approximately the same condition of 
trimmed level flight. The flight condition was between approach and cruise air- 
speed (approximately 40 percent power). For the longitudinal tests, the mean 
Initial trimmed airspeed was 45.3 m/sec with a standard deviation of 1.7 m/sec 
and the initial altitude was 641 m with a standard deviation of 57 m. For the 
lateral tests, the mean initial trimmed airspeed was 45.8 m/sec with a standard 
deviation of 1.6 m/sec, and the initial altitude was 564 m with a standard devi- 
ation of 24 m. Effects of variations in airspeed and altitude were considered 
to be negligible and were Ignored in the processing and analysis of the data. 

It was necessary to apply corrections to some of the data before processing 
with the identification algorithm. Airspeed was corrected for position error 
and altitude (assuming standard temperature) to obtain true airspeed. Then air- 
speed, angle of attack, and angle of sideslip were corrected for upwash and 
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resolved into velocity components u, v, and w by using the relationships 
developed in appendix B. The control-surface positions were measured at the 
respective cables in the vicinity of the cockpit. Because some of the data were 
commutated, it was necessary to reconstruct the data between data points by lin- 
ear interpolation and then sample simultaneously to avoid time-shift errors due 
to commutation. The continuous data were sampled digitally by using a zero- 
phase Ormsby filter (ref. 10). 

The accelerometers were considered to be located on the airplane center of 
gravity in the X-Y plane but their z-position relative to the 'airplane center 
of gravity was unknown. The computer program used for identification had the 
option of identifying the z-position of the Y-accelerometer ; so this variable 
was activated in the algorithm for the appropriate maneuvers. Therefore, the 
accelerometer data were not corrected to the center-of-gravity position before 
being examined by use of the identification algorithm. A brief description of 
the identification algorithm is contained in appendix C. 


CONTROL INPUTS 

Various control inputs were used to determine dependence of parameter con- 
sistency on the inputs. Some of the inputs were suggested by several input 
design studies (refs. 6, 7, and 8) and others were inputs common to flight test- 
ing. The inputs used in this study were not the optimal inputs of these studies 
but were simplified forms. 

For the identification of the longitudinal parameters, three basic types 
of inputs were attempted: the stabilator square wave, sine wave, and rapid 

rise followed by slower decay. The square-wave input (input A, fig. 3) was 
chosen because it was thought to contain the frequency content necessary to 
excite the short-period mode. The period and amplitude (stabilator travel) 
were chosen for ease of pilot implementation, as well as to keep pitch-attitude 
changes within 5° to 10^ of trim. The stabilator sine-wave input (input B, 
fig. 3) is often used to characterize a second-order system. The stabilator 
rapid rise followed by slower decay (input C, fig. 3) was an attempt to approxi- 
mate an input form suggested by reference 6. 

For lateral identification, inputs consisted of rudder or ailerons applied 
individually or a sequential combination of both. (See fig. 4.) Switching-type 
(square-wave) inputs were recommended in several references (for example, refs. 7 
and 8) as an approximation to the optimal input; therefore, square-wave inputs 
were attempted, although the switching times and amplitudes used for the tests 
were not optimal. Both sine- and square-wave forms were investigated for each 
control. The resulting inputs were rudder square wave (input D), rudder sine 
wave (input E), aileron square wave (input F), and aileron sine wave (input G) . 

Rudder inputs alone do not provide adequate excitation of the lateral modes 
for parameter identification (ref. 11). To provide better excitation than single 
controls (rudder or ailerons individually) can produce, combinations of rudder 
and ailerons were used sequentially. That is, the aileron inputs were imple- 
mented followed immediately by rudder inputs or rudder inputs were followed by 
aileron inputs. The resulting sequential inputs were rudder square wave fol- 


7 



lowed by aileron square wave (input H), rudder sine wave followed by aileron 
sine wave (input I), aileron square wave followed by rudder square wave 
(input J), and aileron sine wave followed by rudder sine wave (input K). 

Typical time histories of these inputs can be seen in figure 4. 

Several repeat runs were made for each input in an effort to obtain the 
desired input form. Inspection of the data for each input form showed that none 
of the input repeats excelled over the others, so all the runs were used for 
analysis. For the longitudinal inputs, this resulted in five repeats of input A 
and two of B and C. For the lateral inputs, D and E each had four repeats, F 
had two, G had four, H had five, I had three, J had four, and K had two. 


CRITERIA USED FOR EVALUATING INPUTS 

The results from each of the several different control inputs were examined 
before presenting a set of stability and control derivative values which best 
described the subject airplane. In determining which inputs provided the most 
consistent identification of the parameter values, three criteria were consid- 
ered: ensemble variance of the parameter estimates, estimated Cramer-Rao lower 

bound (ref. 9), and parameter correlations. 

The ensemble variance is the variance of the parameter estimates from the 
same type of input. These variances of the parameter estimates were compared 
by using the F-tests (ref. 12). This comparison indicated which inputs produced 
the most consistent identification of the unknown parameters. 

The next criterion, the estimated Cramer-Rao lower bound, was examined also 
using the F-test for further indication of the effectiveness of the inputs for 
identification. According to reference 9, the estimated Cramer-Rao lower bound 
is proportional to the uncertainty level in the estimate of an unknown parameter. 
The primary function of these uncertainty levels used here was to indicate which 
inputs provided the greatest information and, inversely, the least uncertainty 
in defining a particular derivative value. The higher the uncertainty level for 
a particular derivative, the less information was generated by the input to 
define that derivative. Therefore, a comparison of the uncertainty levels 
obtained for a parameter from various Inputs gave some indication of which 
inputs provided the best identification of that parameter. 

Another criterion that aids in discerning which inputs provide the more 
reliable identification of the unknown parameters is the correlations, or depen- 
dencies, between pairs of parameters (ref. 6). These correlations between param- 
eters are determined from the off-diagonal terms of the error covariance matrix 
for estimated parameters (parameter covariance matrix). Correlation values near 
±1 indicate parameter dependencies, that is, changes in one parameter value caus- 
ing changes in another. These parameter dependencies do not necessarily arise 
from physical relationships, but may also arise from inadequate excitation of 
some of the states for which parameters are to be identified. Inadequate exci- 
tation can be analogous to the problem of attempting to identify too many param- 
eters (ref. 6). High correlations can result in a nonuniqueness of the identifi- 
cation problem which can result in parameter ratios or sums of parameters being 
identified rather than the specific parameters themselves. Although reliable 


8 



parameter estimates can occur even when correlations are high, in many cases non- 
realistic parameter values and significant run-to-run scatter have occurred when 
correlation coefficients are 0.9 or larger. Therefore, large correlation coeffi- 
cients are an indication of potential identification problems. One way these 
parameter dependencies can be reduced is by choosing inputs that adequately excite 
the modes for which the associated parameters are to be identified. Therefore, 
examination of the parameter covariance matrix associated with each input 
revealed which input reduced parameter dependencies or correlations and, in 
turn, probably provided the more reliable parameter estimates. 


RESULTS AND DISCUSSION 

In order to evaluate the uniformity of input repeats, the mean input period 
and amplitude and the associated standard deviations for each input are pre- 
sented in table III. As can be noted, there was considerable variability in the 
mean period among the lateral inputs (5.4 to 10.4 sec). Generally, this varia- 
bility did not appear to correspond to the scatter in the parameter estimates. 


Effects of Various Longitudinal Control Inputs 

The longitudinal parameter values obtained from applying the maximum like- 
lihood technique to the longitudinal flight test data are presented in figure 5 
and table IV. In the upper half of figure 5 are the longitudinal parameter val- 
ues obtained for each run for inputs A, B, and C. Note that the run numbers 
are located beside the extracted values and, in some cases, the extracted values 
overlap so closely that one symbol may represent two or three data points. In 
the lower half of figure 5 are the mean values for each type of input, with the 
upper and lower bars representing plus and minus values of the standard deviation 
for each parameter. Also, just to the right of these bars are bars representing 
the mean plus or minus the mean of the estimated Cramer-Rao lower bound for each 
input. 

Comparing the variances (F-test) of the longitudinal parameter estimates 
across the inputs indicated no significant differences. Therefore, based on 
the variance criterion, none of the stabilator inputs showed significant improve- 
ment in providing consistent estimates of the parameters. 

The next criterion to be considered was the uncertainty level (based on the 
estimated Cramer-Rao lower bound) in the parameter estimates. A visual inspec- 
tion of figure 5 shows little difference in longitudinal parameter uncertainty 
levels among the inputs. The F-test was applied to the square of the uncer- 
tainty levels of the parameter estimates for the three inputs. Based on uncer- 
tainty levels, none of the longitudinal inputs stood out as significantly improv- 
ing the consistency of the parameter estimates. This same conclusion was reached 
by considering the variance, so the conclusions for these two criteria are 
consistent. 

The third criterion, the parameter correlations, was considered next. For 
the present study, basically the only longitudinal parameters which appeared to 
be correlated were Cm and input A, the correlation between these 
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two parameters rainged from 0.90 to 0.95; for input B, from 0.93 to 0.96; and for 
input C, from 0.84 to 0.91. Also, for input B, run 6 exhibited relatively high 
correlations for some of the other parameters also. Based on parameter correla- 
tions, input C appeared to offer slightly improved identification. The form of 
input C was loosely patterned after the optimal input of reference 6, which was 
based on criteria related to the trace and determinant of the Fisher information 
matrix. This matrix is the inverse of the parameter covariance matrix. 

Based on the three criteria considered, none of the three longitudinal 
inputs clearly offered better identification than the other two. 


Effects of Various Lateral Control Inputs 

The lateral parameter values obtained from the flight data are shown in 
figure 6 and table V. The lateral data are presented in the same format as that 
used for the longitudinal data. Note that no attempt was made to identify aile- 
ron control derivatives with rudder inputs or rudder control derivatives with 
aileron inputs. 

A visual inspection of the data reveals that inputs with rudder alone (D 
and E) generally gave the least consistent estimates of the parameters except 
for the rudder control derivatives. Correspondingly, ailerons alone (inputs F 
and G) generally provided more consistent estimates than the rudder alone. For 
one of the aileron runs (run 22 of input G), however, the parameter estimates 
were considerably out of line for the rolling- and yawing-moment parameters. 

An examination of the time history of run 22 showed that the period of the input 
was about 2 sec longer than the other runs for input G, although the basic form 
was the same. This long input period may not have properly excited the airplane 
dynamics and thus may have resulted in some identification problems. Parameter 
correlations for run 22 were exceptionally high and are discussed subsequently. 
Due to the problems encountered with run 22, the associated parameter values 
were omitted from the computation of the statistics shown in the lower half of 
figure 6. 

To investigate lateral inputs as they relate to consistent identification, 
the three criteria - ensemble variance, uncertainty level, and parameter correla- 
tions - were again used as for the longitudinal parameters. 

In order to determine whether input form affected consistency, results 
from inputs D and E, and F and G were compared by using the variance criterion 
(F-test). For rudder input forms, D generally provided more consistent esti- 
mates than E but only Cyg, ^n^ were significantly more consistent 

(80-percent level). In no case did input E provide significantly greater consis- 
tency in the parameter estimates than input D. This result seems to indicate 
that, for rudder inputs, rudder square waves provided the better identification. 
For aileron inputs, F provided more consistent estimates than G (sine waves). 

Input F demonstrates significantly greater consistency for Cvo , Cv , C, ^ , 

p P p 

and Cn^a than input G. Therefore, it appears that square-wave inputs provided 


10 



the more consistent estimates for ailerons also, but uncertainty levels (esti- 
mated Cramer-Rao lower bound) and parameter correlations should be examined 
before a conclusion is drawn. 

Rudder inputs were next compared to aileron inputs. Comparing rudder 
square- wave inputs (D) with aileron square-wave inputs (F) showed that for 
Cyg, Cy^, Cy^, and C-Lg, ailerons provided significantly (80-percent level). 

more consistent identification and in no case did input D provide more consis- 
tent estimates than any of the runs using aileron inputs. A similar trend, 
although less pronounced, held when comparing inputs E and G (rudder and aile- 
ron sine waves), with Cy and being significantly more consistently 

p p 

identified for the aileron inputs. It is interesting to note that for single 
control inputs, the derivatives which might be assumed to be associated with 
excitations from the rudder inputs appear to be better defined by using aileron 
inputs. One explanation is that the correlations between parameters were higher 
for the rudder inputs, making identification of the correlated parameters diffi- 
cult. Correlations are discussed in more detail subsequently. 

Upon comparing the single control inputs (D, E, F, and G) with the sequen- 
tial inputs (H, I, J, and K) it was found that the sequential inputs provided 
significantly more consistent estimates of most of the stability derivatives, 
with the exception of input F. Input F was nearly as consistent as the sequen- 
tial inputs, but this may be due to the small data sample (2 runs). For the 
control derivatives, the single inputs generally provided estimates which were 
equally as consistent as those obtained with sequential inputs. 

Inputs H and J, and I and K were compared to determine whether the order 
of the inputs for the sequential inputs was significant. Input J gave signifi- 
cantly (80-percent level) more consistent estimates of the rolling-moment deriva- 
tives except for input H gave significantly more consistent estimates 

of Cy^ and These trends did not hold when comparing inputs I and K; 

only identified as significantly more consistent and this was for 

input I. Therefore, based on the information examined thus far, it is difficult 
to conclude whether order makes any difference. 

The next criterion to be considered in evaluating input effectiveness for 
identification purposes was the uncertainty level, or estimated Cramer-Rao lower 
bound. Generally, uncertainty levels were less than the standard deviations, 
especially for the single inputs (D, E, F, and G). Greater uncertainty levels 
in many cases corresponded to greater standard deviations, which indicated that 
uncertainty levels were somewhat analogous to standard deviations, as stated in 
reference 9. Therefore, the trends determined by using the F-test on the vari- 
ances should generally hold when considering the square of the uncertainty 
levels. Looking at input form (comparing data from inputs D and E, and F and G) 
the distinction between parameter estimates that was noted when comparing ensem- 
ble variance was not as obvious when comparing the square of the uncertainty 
level. Statistically speaking, none of the estimated parameters had signifi- 
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and with coefficient values ranging from 0.85 to 0.89 and C^p and 

with values ranging from 0.815 to 0.869. Inputs I and K were not quite as effec- 
tive in reducing correlations as H and J but were considerably more effective 
than D, E, F, and G. 

Based on the information contained in the parameter covariance matrix, the 
sequential inputs clearly reduce parameter correlations, or dependencies, compared 
to the single control inputs. This implies that the more reliable lateral iden- 
tification is obtained from the sequential inputs. 

The z-position of the Y-accelerometer was also identified for the lateral 
runs. As for the parameter estimates, the most consistent estimate of the loca- 
tion was observed for the sequential inputs. These estimates ranged from 0.24 
to 0.30 m below the airplane center of gravity. 


Use of Estimated Parameters in Mathematical Model 

None of the longitudinal inputs attempted offered clear improvement in the 
consistency of the extracted derivatives. Therefore, the resulting parameters 
presented for the subject airplane are the arithmetic means of the various 
parameter values extracted for all 9 longitudinal runs. These mean values 
along with standard deviations are presented in table VI. 

For the lateral data, the sequential inputs generally provided the most 
consistent estimates of the lateral parameters. Therefore, the arithmetic means 
of the parameters obtained from the sequential inputs were used to describe the 
airplane in the subsequent analysis. These mean values are also presented in 
table VI. 

The mean parameter values of table VI were used to generate the computed 
time histories of figures 7 and 8. The match between the computed airplane 
responses and the measured responses was considered to be good, although a 
slightly better fit could be obtained with the actual values for that run. 

The last column in table VI contains the standard deviations as a percent 
of the mean value for each derivative. This quantity, the coefficient of varia- 
tion, provides a measure of the relative consistency of one derivative with 
another. For example, (3.6 percent) was the most consistent lateral deriv- 

p 

ative identified and (4.0 percent) was the most consistent longitudinal 

derivative. Generally, the static derivatives such as and were more 

consistently identified than were the rotational and control derivatives. 


Comparison of the Extracted Parameters 


One question which arises whenever derivatives are extracted from flight 
data is. Are the results reasonable? Good agreement between measured and com- 
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puted time histories does not necessarily guarantee that the parameters are 
reasonable. Correlations, an inadequate mathematical model, or other problems 
could exist which could drive the parameters to unrealistic values to fit the 
data. Therefore, to lend some confidence to the extracted parameter values, 
these values were compared with results for the same airplane, which were deter- 
mined by another method. The second column of table VII contains the estimated 
parameter values given in table VI. The third column of table VII contains an 
estimate of the stability and control derivatives of the subject airplane 
obtained by using the theoretical-empirical techniques presented in refer- 
ence 13. These values were transformed to body-axis derivatives by the trans- 
formations presented in reference 14. Upon comparison, it can be seen that the 
estimated parameter values of the present study are generally consistent with 
those computed from reference 13. 

As a matter of interest, the derivative values obtained by using the maxi- 
mum likelihood technique, as well as those from reference 13, were used in the 
equations of motion to generate time histories. The computed time histories 
based on the derivative values and the accelerometer offset obtained by using the 
maximum likelihood technique and the measured flight data time histories are 
shown in figures 7 and 8. The computed time histories based on the derivatives 
obtained by using the methods of reference 13 and the same accelerometer offset 
and the same flight data time histories are shown in figures 9 and 10. As would 
be expected, the maximum likelihood derivatives provided the better fit to the 
flight data, which was especially noticeable for the lateral data. 


CONCLUDING REMARKS 

Stability and control derivative values were determined by use of the maxi- 
mum likelihood algorithm from flight test data taken in cruise conditions for a 
low-wing, single-engine airplane. Various control inputs were used to determine 
whether certain inputs improved the run-to-run consistency of the derivative 
estimates. The longitudinal inputs included stabilator square waves and sine 
waves and a rapid rise followed by a slower decay. The lateral inputs consisted 
of rudder and aileron inputs (square and sine waves) individually and in 
sequence. 

The parameter values obtained from the various inputs were compared for 
improved consistency by using three criteria: the ensemble variance of the 

parameter estimates, the estimated Cramer-Rao lower bound, and parameter corre- 
lations. The results indicated that none of the longitudinal inputs resulted 
in improved identification. On the other hand, the lateral inputs which con- 
sisted of rudder and ailerons in sequence clearly offered improved consistency 
over the rudder or aileron inputs individually. 

For the single inputs, the square-wave inputs for both rudder and ailerons 
generally provided the more consistent parameter estimates when compared with 
their sine-wave counterparts. For the sequential inputs, square-wave inputs 
appeared to provide more precise estimates of the parameters. The order of the 
sequential inputs made little difference in the consistency of the estimates. 
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In order to validate the parameter values determined by use of the maximum 
likelihood algorithm, parameter values for the subject airplane were also deter- 
mined by a theoretical-empirical method. These values and those obtained by 
using the maximum likelihood algorithm were compared both individually and in 
their ability to describe the airplane dynamics accurately. The results indi- 
cated that the values from the maximum likelihood algorithm more nearly described 
the airplane motions. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 17, 1 977 
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APPENDIX A 


EQUATIONS OF MOTION 


The equations used in this program are perturbation equations from trimmed 
level flight and are written relative to the set of body axes shown in figure 1. 


The equations used to describe the longitudinal motions were 
u = -qw + rv - g sin 0 + 1 p “ *^t^J 

w = -pv + qu + g cos 0 cos 4> + I p C 2 © “ *^t^ ^Zn ^ 

2 m |_ ^ H ?v 


(AD 


q = pr + ^(r2 _ p2) + p + Ci^(a - at) + |£ 

ly ly 2Iy L ^ 


(A2) 


+ Cnjq " ^e,t) 


0 = q cos 4> - r sin (j) 


ay = l(u + qw - rv + g sin 0) 
S 


(A3) 

(A4) 

(A5) 


a-Z = — (w + pv - qu - g cos 0 cos 4>) 


(A6) 


V = \/u2 + 


(A7) 

a = tan"^ 

w 

(A8) 


u 


a ~ 


(A9) 

u 



^Z;c “ — 

^®6e 

(A10) 
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APPENDIX A 


The values of the lateral states v, p, r, and ^ used in the longitudinal 
equations were the flight-measured quantities. 

Since thrust changes are not explicitly modeled in the equations of motion, 
^Xrw equation (A1) and of equation (A2) are not necessarily pure Cx^ 

vJC Ob Ob 

and but may contain small contributions due to changes in thrust. There- 

fore, Cx^ and ^ 2 ^^, as determined in this study, are given by 

Cy^ = cos e 

aa 


07^ = + C*p sin e 

Ta 


Since, in this study, thrust was held constant and the angle-of-attack changes 
were no more than 7^ peak to peak, the contributions of thrust to Cx and 
were considered minimal. 


The equations used to compute the lateral motions were 

p ^ 

2 m L 


V = -ru 


+ pw + g cos 9 sin <j> + 1 P o + Cy — 

o m I ^ p r 2V 


+ Cv 

pb 

2V 

+ ( 

Ixz 

+ 1 

^ly 

Ix 

\ 

( ■ 

+ c * 1 ^ 

rb 

+ ( 

Lr 

2V 


Ixz p 

+ j 

fix 

Iz 

1 

1 ^ 

+ 

rb 

+ ( 

xir 

2V 



(All) 


ix } 


Ix 


' 1,0 




P 2V 


'•6a 


(A12) 


pq 


- (M)Qr . 1 p . C„ 


Pb 

P 2V 


<j) = p + (q sin 4) + r cos (j)) tan 9 


ay = 2(v + ru - pw - g cos 9 sin 4>) 
g 


(A13) 

(A14) 

(A15) 
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V = + v^ + w^ 

(A16) 

3 = sin-1 V 

(A17) 


V 


The values of longitudinal states u, w, q, and 0 used in the lateral equa- 
tions were the flight-measured quantities. The equations were used to compute 
the airplane state responses. The computed responses were then compared with 
the recorded responses from the flight tests and the differences were used to 
update the parameters (stability and control derivatives) to improve the fit. 

The longitudinal measured and computed responses, or states, used in the 
algorithm for this study were u, w, q, 0, ax, and a2- The lateral states 
used were v, p, r, (j), and ay- An abbreviated discussion of the identifica- 
tion algorithm is given in appendix C. 



APPENDIX B 


TRANSFORMATION OF MEASURED V, a, AND B 

The boom on which the dynamic pressure, angle of attack, and angle of side- 
slip were measured was located at the left wing tip parallel to the airplane 
X body axis. The sensing elements were located about 3/4 c ahead of the lead- 
ing edge of the wing. These measurements were corrected for upwash and trans- 
formed to the airplane center-of-gravity position. The following development 
describes the procedure used for the data of this study. The following sketches 
show the side view 



and the top view 


Y 



of the schematic of the left wing tip. The equations related to these sketches 
are 


Uffl = 

Vm 

cos 



Wm = 

Vm 

sin 



% = 

Urn 

tan 


(B1) 

Vm = 

Urn 

tan 


(B2) 


The magnitude of the velocity vector was measured by the pitot static tube and 
the direction was given by the a and 3 vanes. The magnitude of the velocity 
vector is given by the equation 
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Vm = \/%^ + Vn,2 + (B3) 

Substituting equations (B1) and (B2) into equation (B3) gives 

(B4) 

+ tan^ 3iii + tan^ Om 

It was assumed that the v-component was unaffected by the upwash and needed 
no upwash correction. Therefore, the problem was considered to be contained in 
the X-Z plane as shown in the following sketch: 




As can be seen from this sketch, in order to correct velocity components Um 
and WjQ for upwash, the affected velocity vector needs to be rotated an amount 
which will result in u' and w* being the correct components when the rotated 
vector is projected onto the X- and Z-axes. 

Based on experimental data, the angle of attack corrected for upwash a 
was determined to be 0.75<%; that is, a = 0.750^. Therefore, 

•J' = cos (0.75Onj) 

V’ = Vm 


w’ = \J^m^ + sin (0.75ajn) 

Substitution of equations (B1), (B2), and (B4) into these equations yields 

ui ^ VmVl + tan^ (% cos (0 .75am) 
y1^Ttan^Bm + tan^ 

= Vm tan 3m 

y1 + tan^ Bjjj + tan^ am 
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_ VmVl + tan^ am sin (0.75am) 

\j'l + tan^ gjjj + tan^ am 

The components u' , v', and w' are the upwash corrected components of 

velocity at the boom. These components must also be corrected to the airplane 
center-of-gravity position by removing velocities resulting from rotational 
rates about the center of gravity. This correction results in the following 
equations : 


u = Ujn cos (O.TSotm) - qz + ry (B5) 

V = Um tan 3m “ 

w B Um sin (0.75am) - py + qx (B7) 

where x, y, and z denote the position of the sensors relative to the air- 
plane center of gravity. 
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APPENDIX C 


PARAMETER ESTIMATION PROCEDURE 

The parameter estimation procedure used in this study is knovm as the maxi- 
mum likelihood technique. This formulation of the maximum likelihood technique 
considered only measurement noise. The technique utilizes the differences 
between the measured and computed states (fit error) to estimate the unknown 
values of the parameters necessary to minimize the fit error. The measurement 
noise 


ri(ti) = xM(ti) - Xc(a,ti) 

is assumed to be white with Gaussian distribution. This assumption is used to 
construct the likelihood function, the natural logarithm of which is 

N 

In L(aO,Ri) = -I ^ fiT(ti)Ri-'' q(ti) - N In |Ri 1 
^ i=1 ^ 


where R-| is the estimate of the state error covariance matrix based on the 
nominal solution, t^ is time, and otO is the nominal parameter vector. 

To estimate the parameter vector a the likelihood function was maximized 
with respect to 3° and Ri. Because of the nonlinear nature of the estimation 
problem, it was necessary to use an iterative procedure to estimate the parame- 
ters. Therefore, x(a,ti), the estimated state vector at time tj^, vras expanded 
about the nominal parameter vector, neglecting second-order and higher order 
terms. That is. 


x(5° + Aa,ti) = x(a°,ti 



o 

A«j 


expresses the estimated state vector as a function of the nominal solution plus 
the change in the state vector estimate as a function of the computed change in 
the parameter vector. This expression for the state vector estimation is then 
used in the fit error vector to obtain 


n(ti) = 5 ^M(ti) - x°(ti) - A(ti)°Aa = n(ti)° - A(ti)°Aa 


where 



and defines state sensitivity to the current parameters. Therefore, the natural 
logarithm of the likelihood function becomes 
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In L(3 o + Aa,R-|) 


N f 

-1 ^ - A(ti)Ao]''' R-i-IpCti) 


- A(ti)Aa 



Maximizing the likelihood function with respect to Aa, 
3 In L(gO + A^,Ri) _ q 
3 ^ 


I in |Ril 


and solving for Aa yield the parameter change equations, 



N 

-1 

' N 

Aa = 

XI AT(ti)Ri-"' A(ti) 


AT(ti)Rl-'’ n(ti) 


_i=1 


_i=1 


which are used to compute the Aa that tends to maximize the likelihood func- 


tion. Note that the matrix 


N 


^ AT(ti)Rr^ A(ti) 


i = 1 


-1 


is the estimate of the 


error covariance matrix for the estimated parameters (parameter covariance 
matrix). This matrix contains the parameter dependencies or correlations which 
are used to construct parameter correlation coefficients - 


To estimate R^, the measurement noise covariance matrix, the likelihood 
function can be maximized with respect to R-j , which is approximately 

N 

Ri = I XI 

according to reference 1 • 

The effect of maximizing the likelihood function is the same as minimizing 
the cost function 


r ^ 

j = det (1 n(ti)n'^(ti 

In using the program, the cost function was computed and displayed on the pro- 
gram operator console (ref. 4). Generally, with each iteration the fit improved, 
as evidenced by a reduction in the cost function. Once the cost function settled 
out so that changes from iteration to iteration, as defined by (Jj^ - Jk+i)/Jk, 
were less than 0.01, the parameters which maximize the likelihood function or 
minimize the cost function (fit error) were considered identified. The expres- 
sion (Jk - Jk+1^/*^k represents the change in the fit error from successive 
iterations divided by the fit error of the previous iteration. To avoid' possi- 
ble correlation problems (see section "Criteria Used for Evaluating Inputs") 
between and 0^^ was held fixed at -4.00 (estimated in a preliminary 

study by using ref. 13), so essentially the program was identifying the combina- 
tion Cjjq + Cnj^ with Cn,^ held constant. 
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Also, 


to further reduce correlation problems between estimated parameters, 
Cmgg were assumed to be geometrically related as 


^Z6e = ^ ^“6e 


as was done in reference 4 . The program described in reference 1 was modified 
so that was not treated as an active parameter. Then Cnj.. was 


extracted, and ^2^^ was calculated for each iteration and its value was sub- 


stituted into the equations of motion. Although ^2^^ was not active, its 


value was changing from iteration to iteration, 
in sensitivity equations, £_ 

t 


To account for this variation 
was substituted for ^2^^ in 


the w equations and then the partial derivatives were taken of the w equa- 
tion with respect to and ^2gg. 
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TABLE I.- GEOMETRIC DETAILS OF SUBJECT AIRPLANE 


Mass, kg . . 1074.1 

Inertia: 

l x, kg-m2 1220 

ly, kg-m2 1898 

ly, kg-m2 2712 

IXZ> kg-m2 68 

Fuselage length, m 6.9 

Wing: 

Area , 14.9 

Aspect ratio 5.625 

Span ,m 9.1 

Mean geometric chord, m 1.62 

Vertical tail: 

Area, 1,07 

Aspect ratio 2.02 

Span, m 1.47 

Rudder area, 0.38 

Horizontal tail (stabilator) ; 

Area, m^ 2.32 

Aspect ratio 4.0 

Span, m 3. 05 

Tail length, center of gravity to quarter-chord point of mean 

geometric chord, m 4.2 

Location relative to airplane center of gravity of sensors to measure 
velocity, angle of attack, and angle of sideslip; 

X, m 1.29 

y, m _4.46 

z, m.................. Q 
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TABLE II.- INSTRUMENTATION RANGES 


Instrument function 

Range 

Airspeed, m/sec 

0 to 61 

Angle of attack, deg 

-30 to 100 

Angle of sideslip, deg 

-60 to 60 

Altitude, m 

0 to 3050 

Normal acceleration, g units 

-3 to 6 

Longitudinal acceleration, g units 

-1 to 1 

Lateral acceleration, g units 

-1 to 1 

Elevator position, deg 

Full throw (1.5 to -19-0) 

Aileron position, deg 

Full throw (-43 to 43) 

Rudder position, deg 

Full throw (-25 to 25) 

Throttle position 

Full throw 

Engine speed, rpm 

0 to 2700 

Pitch rate, deg/sec 

-100 to 100 

Roll rate, deg/sec 

-180 to 180 

Yaw rate , deg/sec 

-180 to 180 

Pitch attitude, deg 

-85 to 85 

Roll attitude, deg 

-180 to 180 
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TABLE III.- CONSISTENCY OF INPUTS USED IN FLIGHT TEST 


Input 

Control 

Input form 
(a) 

Mean 

period, 

sec 

Standard 
deviation 
of mean 
period, 
sec 

Mean peak-to- 
peak amplitude, 
deg 

Standard 
deviation 
of mean 
amplitude , 
deg 

A 

Stabilator 

Square wave 

4.2 

0.6 

2.4 

0.4 

B 

Stabilator 

Sine wave 

5.5 

1.0 

2.5 

.4 

C 

Stabilator 

Rapid rise, slow 
decline each 
direction 

4.4 

.8 

4.8 

.4 

D 

Rudder 

Square wave 

5.5 

.7 

11.2 

1.0 

E 

Rudder 

Sine wave 

6.6 

1.2 

20.8 

2.8 

F 

Ailerons 

Square wave 

5.9 

.6 

1^15.5 

2.1 

G 

Ailerons 

Sine wave 

8.0 

1.1 

bl5.8 

1.5 

1 

H 

Rudder and 

1 

Rudder square 

5.4 

.5 

10.0 

1 .0 


ailerons 

wave followed 
by aileron 
square wave 

5.9 

.6 

bl7.1 

2.2 

I 

Rudder and 

Rudder sine wave 

7.9 

.8 

15.7 

2.9 


ailerons 

followed by 
aileron sine 
wave 

8.1 

1.2 

t^is.o 

2.6 

J 

Ailerons and 

Aileron square 

5.8 

.5 

>^15.3 

1.7 


rudder 

wave followed 
by rudder 
square wave 

6.0 

.4 

12.4 

1.5 

K 

Ailerons and 

Aileron sine wave 

8.8 

1.1 

t>15.1 

1.3 


rudder 

followed by 
rudder sine 
wave 

10.4 

.8 

17.5 

3.5 


^Inputs were not strictly sine waves or square waves but were pilot 
attempts at these forms. See figures 3 and 4 for typical input time histories. 
“Aileron deflections are defined as left aileron minus right aileron. 
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TABLE IV.- ENSEMBLE MEAN AND STANDARD DEVIATION FOR EACH LONGITUDINAL 
INPUT ALONG WITH ESTIMATED CRAMER-RAO LOWER BOUND 


Parameter 

Mean value 

Standard deviation 
of mean value 

c^CRLB 

®'^CRLB 



Input A (5 runs) 



Cx,o 

0.0533 

0.0068 

0.00022 

0.00010 

Cx 

Aa 

.6481 

.2275 

.01830 

.00735 

Cz,o 

-.5533 

.0465 

.00060 

.00029 

C7 

ZlQt 

-4.3828 

.1864 

.03826 

.00947 


- 15.1932 

4.1709 

.80940 

.22418 


-.4394 

.0324 

.00356 

.00148 

0 

B 

-8.7746 

3.7690 

.20806 

.09952 


-1.5265 

.2125 

.01830 

.00870 


Input B (2 runs) 



^X,o 

0.0405 

0.0156 

0.00013 

0.00007 

Cx 

.7424 

.1573 

.00815 

.00233 

Cz,o 

-.5829 

.0573 

.00046 

.00024 


-4.2332 

.0433 

.02680 

.00424 


-20.9335 

1.0257 

.59945 

.06074 


-.4723 

.0071 

.00260 

.00057 

*^mq 

-7.6263 

2.6367 

.16445 

.01407 


-1.5265 

.2850 

.01560 

.00042 
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TABLE IV.- Concluded 


Parameter 

Mean value 

Standard deviation 

‘^CRLB 

® Op pin 



of mean value 


OriLo 



Input C (2 runs) 



Cx,o 

0.0^69 

0.0052 

0.00020 

0.00006 

CXa 

.5076 

.1623 

.02120 

.01018 

Cz,o 

-.5644 

.0508 

.00054 

.00009 

^Za 

-4.4514 

.2157 

.04105 

.00615 


-17.0211 

5.2584 

.72725 

.04632 


-.4925 

.0311 

.00390 

.00028 

^mq 

-8.4654 

1 . 1422 

.16715 

.04674 

^®6e 

-1.5786 

.0947 

.01500 

.00311 



TABLE V.- ENSEMBLE MEAN AND STANDARD DEVIATION FOR EACH LATERAL INPUT 


ALONG WITH ESTIMATED CRAMER-RAO LOWER BOUND* 
(a) Input D (4 runs) 


Parameter 

Mean value 

Standard deviation 
of mean value 

^CRLB 

^‘^CRLB 


-0.5737 

0.0799 

0.0372 

0.0229 

% 

-.0291 

.5170 

.2723 

.1737 

Cy^ 

.2806 

.2652 

.0730 

.0105 

%r 

.0342 

.0090 

.0040 

.0007 

^^3 

-.0631 

.0229 

.0030 

.0021 


-.3400 

.1375 

.0179 

.0127 

"ir 

.0376 

.0475 

.0060 

.0034 


.0029 

.0038 

.0006 

.0002 

^n3 

.0521 

.0098 

.0018 

.0008 

Cnp 

- . 0480 

.0619 

.0107 

.0045 

Cn, 

-.0923 

.0212 

.0034 

.0012 


-.0338 

.0021 

.0003 

.00004 



TABLE V.- Continued 


(b) Input E (4 runs) 


Parameter 

Mean value 

Standard deviation 

‘^CRLB 

®°CRLB 

0.0100 

Cyg 

-0.4639 

of mean value 
0.1868 

0.0332 

Cv 

.8892 

1.0266 

.1879 

.0625 


.7896 

.2648 

. 1408 

.0083 


.0764 

.0210 

.0074 

.0008 


-.0438 

.0174 

.0024 

.0008 

Cl 

-.2110 

.0907 

.0132 

.0057 

Cl 

.0647 

.0526 

.0076 

.0034 

P, 

.0068 

.0018 

.0007 

.0003 

Cn3 

.0492 

.0203 

.0023 

.0006 

Cnp 

-.0802 

.1216 

.0130 

.0038 

Cn, 

-.1042 

.0502 

.0052 

.0024 

c „ 

-.0368 

.0014 

.0004 

.0002 

^n6a 







TABLE V.- Continued 


(c) Input F (2 runs) 


Parameter 

Mean value 

Standard deviation 
of mean value 

<^CRLB 

®°^crlb 


-0.5396 

0.0001 

0.0061 

0.0008 

CYp 

.0947 

.0079 

.0162 

.0012 

CYp 

.3182 

.0339 

.0265 

.0031 






-.0413 

.0010 

.0012 

.0003 


-.2538 

.0262 

.0081 

.0001 

"ir 

c 

.0479 

.0245 

.0043 

.0018 


-.0390 

.0018 

.0009 

.0001 

Cng 

.0598 

.0038 

.0004 

.0008 

Cn 

np 

-.0462 

.0101 

.0037 

.0014 

Cn 

“r 

c 

-.1072 

.0145 

.0022 

.0004 

^n6a 

.0058 

.0002 

.0004 

.0001 



TABLE V.- Continued 


(d) Input G (3 runs) 


Parameter 

Mean value 

Standard deviation 
of mean value 

‘^CRLB 

^c^CRLB 

% 

-0.5022 

0.0882 

0.0177 

0.0070 

% 

.0328 

. 1789 

.0432 

.0178 

CYp 

. 1250 

.1167 

.0500 

.0137 






-.0284 

.0147 

.0020 

.0010 

Cl 

-.1837 

.0468 

.0101 

.0013 

Cl 

^r 

.0071 

.0473 

.0046 

.0020 

^5r 






-.0288 

.0092 

.0011 

.0002 

Cn3 

.0606 

.0011 

.0009 

.0006 

Cr, 

np 

-.0232 

.0107 

.0064 

.0033 

Cn 

Hr 

-.1008 

.0169 

.0032 

.0008 

^6r 






.0083 

.0022 

.0007 

.0002 



TABLE V.- Continued 


(e) Input H (5 runs) 


Parameter 

Mean value 

Standard deviation 
of mean value 

^CRLB 

®<^CRLB 


-0.5442 

0.0109 

0.0048 

0.0009 


.1132 

.0348 

.0152 

.0030 

Cyr 

.3515 

.0315 

.0280 

.0045 

^Y6r 

.0307 

.0022 

.0042 

.0007 


-.0475 

.0016 

.0012 

.0002 


-.2374 

.0128 

.0056 

.0011 

"tr 

.0766 

.0315 

.0042 

.0007 

^^6r 

.0081 

.0031 

.0008 

.0001 

*^i6a 

-.0397 

.0026 

.0008 

.0002 


.0560 

.0013 

.0004 

.0001 

% 

-.0347 

.0046 

.0021 

.0005 

"nr 

-.0938 

.0118 

.0016 

.0003 

"ngr 

-.0338 

.0015 

.0003 

.0001 

"n6a 

.0068 

.0012 

.0068 

i 

.0012 


TABLE V.- Continued 


(f) Input I (3 runs) 


Parameter 

Mean value 

Standard deviation 
of mean value 

'^CRLB 

®'^CRLB 

% 

-0.5706 

0.0233 

0.0076 

0.0007 

% 

.1786 

.0274 

.0242 

.0033 


.3741 

.0902 

.0396 

.0054 


.0566 

.0179 

.0060 

.0006 

^13 

-.0408 

.0040 

.0016 

.0002 

Cl 

-.2245 

.0125 

.0074 

.0020 


.0678 

.0284 

.0052 

.0008 


.0032 

.0046 

.0009 

.0002 


-.0349 

.0011 

.0010 

.0002 

Cng 

.0553 

.0038 

.0008 

.0001 

^np 

-.0737 

.0071 

.0043 

.0006 

Cn, 

-.0911 

.0169 

.0023 

.0003 

*^«6r 

-.0386 

.0028 

.0004 

.0001 

1 

.0018 

,0031 

.0006 

.0001 



TABLE V.- Continued 
(g) Input J (4 runs) 

®‘^CRLB 
0.0004 


Parameter Mean value Standard deviation ctcRLB 

of mean value 

Cvn -0.5516 0.0091 0.0048 




TABLE V.- Concluded 


(h) Input K (2 runs) 


Parameter 

Mean value 

Standard deviation 
of mean value 

‘^CRLB 

®CfCRLB 


-0.5868 

0.0035 

0.0078 

0.0002 


.1146 

.0267 

.0219 

.0012 


.2929 

.0457 

.0405 

.0009 

^Y6r 

.0623 

.0017 

.0058 

.0001 


-.0483 

.0032 

.0015 

.0002 

Cl 

-.2054 

.0302 

.0056 

.0004 


.0520 

.0056 

.0049 

.0007 


.0106 

.0060 

. 0008 

.0001 


-.0364 

.0025 

.0009 

.00002 

Cn3 

.0566 

.0008 

.0006 

.00001 

Cn 

np 

-.0451 

.0113 

.0030 

.0002 

Cn, 

-.0948 

.0035 

.0021 

.0002 


-.0380 

.0013 

.0004 

.00001 


. 0042 

.0008 

.0005 

.00002 



TABLE VI.- MAXIMUM LIKELIHOOD ESTIMATES OF STABILITY AND 
CONTROL DERIVATIVE VALUES FOR SUBJECT AIRPLANE 


Derivative 

Sample mean, 
P 

Sample standard 
deviation, s 

Coefficient of 
variation, s/P x 100 


0.638 

0.198 

31.0 • 

CZa 

- 4.365 

.173 

4,0 

C 7 

-16.875 

4.264 

25.3 

^Z6e 

-.594 

.072 

12.1 

Cma 

-.458 

.035 

7.6 


a-4.000 



''mq 

-8.451 

2.893 

34.2 


-1.538 

. 185 

12,0 

% 

-.558 

.020 

3.6 

Cy 

.124 

.039 

31.4 

CYp 

.370 

.073 

19.6 


.045 

.016 

34.6 


-.046 

.004 

7.9 

Cl 

-.233 

.018 

8.0 


.071 

.023 

32.0 


.006 

.004 

62.3 


-.038 

.003 

7.1 

Cng 

.056 

.002 

3.6 

Cn 

np 

-.048 

.016 

32.8 

Cn, 

-.096 

.011 

11.6 


-.036 

.003 

8.3 


.005 

.002 

40.0 


^Held constant . 



TABLE VII.- COMPARISON OF DERIVATIVES FROM THE PRESENT STUDY 
WITH THOSE FROM METHODS OF REFERENCE 13 



^Held constant. 

^The derivatives this source were mul- 


tiplied by -0.5 to make them compatible with the conventions used 
in present study. 


40 




Figure 1.- System of body axes showing positive sense of velocities, forces, 

and control deflections. 
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Figure 3.- Time histories of stabilator deflections illustrating the three 
longitudinal input forms used in this study. 
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Figure 4,- Concluded. 
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Figure 5.- Parameter values extracted for each longitudinal input along with 
corresponding run numbers and bars depicting ensemble variations and esti- 
mated Cramer-Rao lower bound on variations. 
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Figure 5.- Continued. 


47 


0 r— 


-12 


-16 


-20 


-4 


Run 

06 


-8 

— 


07 

o^ 

08 


0 |— 
-.4 
-.8 


m 


6e 






B 


e/ sq 


'e/si 

Input 


e/ pul 


- 1.2 


- 1.6 


- 2.0 


Run 

2 

5 

3 

1 




e/sq 


-J.- 


Mean ± Standard deviation 

Mean ± Mean of estimated Cramer-Rao lower bound 


06 

07 

_ L _ 

B 


re/ si 
Input 


C 

("^e) pul 


0 





0 





-4 

— 

“T 

1 

“T* 

1 


-.4 

— 




-8 


1 

-4- zjz 

1 

4- =1= 

1 

- 7 - 

*T 

-.8 

C 

— 






1 

1 

-X- 

« -1 

”'6e 





-12 

— 

1 



-1.2 

— 

“T* 

1 









1 


"T” 

-16 





-1.6 


I 

1 

1 

i 

1 


-20 


1 

1 

1 

-2,0 


1 

1 

1 


A 

B 

c 



A 

B 

C 



CD 

Wsl 

('^e) pul 



(6 ) 

\ e/ sq 


e) pul 




Input 





Input 



Figure 5.- Concluded. 
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Figure 6.- Parameter values extracted for each lateral input along with the 
corresponding run numbers and bars depicting ensemble variations and esti- 
mated Cramer-Rao lower bound variations. 
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Figure 6.- Continued 
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Figure 6.- Continued. 
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Figure 6.- Concluded. 
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Figure 7.- Concluded. 
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